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ABSTRACT 

We analyse the covariance of the one-dimensional mass power spectrum along lines of 
sight. The covariance reveals the correlation between different modes of fluctuations in the 
cosmic density field and gives the sample variance error for measurements of the mass power 
spectrum. For Gaussian random fields, the covariance matrix is diagonal. As expected, the 
variance of the measured one-dimensional mass power spectrum is inversely proportional to 
the number of lines of sight that are sampled from each random field. The correlation between 
lines of sight in a single field may alter the covariance. However, lines of sight that are sampled 
far apart are only weakly correlated, so that they can be treated as independent samples. Using 
A^-body simulations, we find that the covariance matrix of the one-dimensional mass power 
spectrum is not diagonal for the cosmic density field due to the non-Gaussianity and that the 
variance is much higher than that of Gaussian random fields. From the covariance, one will be 
able to determine the cosmic variance in the measured one-dimensional mass power spectrum 
as well as to estimate how many lines of sight are needed to achieve a target precision. 

Key words: cosmology: theory - large-scale structure of universe 



1 INTRODUCTION 

The one-dimensional mass power sp ectrum (PS ) and its relation to the th ree-dimensi onal mass PS have been frequently utilized to recover the 
linear mass PS from the Lya forest ICroft et al.l998..1999, 20Q2:.Gnedin & Hamilton.2002. : .,Kim et aU2004i) . This opens a great window 
for studyin g the large-scale structure of the universe oyer a wide range of redshift. As one starts to attempt precision cosmology using the 
Lya forest iCroft et aJ2002l:lMandelbaum et alhoOgLlspergel et all2003l) . it becomes necessary to quantify systematic uncertainties of the 
PS analysis including the covariance. 

For an ensemble of isotropic fiel ds, the one-dimensional mass PS Pio{k) is a simple integral of the three-dimensional mass PS Pw{k) 
iLumsden. Heavens & PeacoclJl98'^ . 



PlD(fc) = 



2-K 



Pao{k')k' dk', 



(1) 



where k is the line-of-sight wavenumber However, there is only one observable universe. One has to replace the ensemble average with a 
spatial average. For instance, one may sample multiple line-of-sight densities from the three-dimensional cosmic density field and use the 
average PS of the one-dimensional densities in place of the ensemble-average quantity. 

The covariance of the spatial-average PS may differ from that of the ensemble-average PS for at least two reasons. First, lines of sight 
(LOS's) are no longer independent of each other. Correlations between close LOS's could inc rease the covariance of the one-dimensional 
mass PS, although they may also be used to extract cosmologic al parameters from the Lya forest iHui. Stebbins & BurlesligigllMcDonald & Miralda-Escudd 
Il999l:lviel et all200ilMcDonalJ2003HRollinde et alhoOgl) . Second, the length of each LOS is always much less than t he size of the uni- 
verse , so that false correlations between different modes are introduced in the covariance (for the three-dimensional case, see iFeld man. Kaiser & Peacockl 
I1994I) . Because the cosmic density field is highly non-Gaussian on scales of interest, A'^-body simulations are helpful for quantifying the co- 
variance. 

The rest of the paper is organized as follows. Section |2| briefly describes the notation and the convention of Fourier transforms, the 
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three-dimensional mass PS, and its covariance. The spatial-average one-dimensional mass PS and its covariance are derived in Sections |3| 
and|4] respectively. The effects of the line-of-sight length is discussed in Section|5| Section|6|presents numerical results of the covariance 
from Ai'-body simulations with different box sizes. The conclusions and discussions are given in Section^ Unless specified otherwise, the 
PS refers to the mass PS. 



2 THREE-DIMENSIONAL POWER SPECTRUM 

We assume the following convention of Fourier transforms for a cubic density field of volume V = B'^: 

l{n)= / 5(a;)e-^"'"-"''-^dx, = \z ^ <5(n) e"™'"/^, (2) 

Jv 

TI — — OC 

where 5(x) and 5(n) are the overdensity and its Fourier counterpart, respectively, the summation y^°° is an abbreviation for y^°° 

I'l./' ^ — /ri — — oo £■ — ^ni ,7i2 ,T13 — - 

and the wavevector fe = 2Tin/B. With the understanding that Fourier modes exist only at discrete wavenumbers, i.e. ni, n2, and 713 are 
integers, one may use k and n interchangeably for convenience. To be complete, the orthonormality relations are 



1 

V ly" — y 



where ^, is the three-dimensional Kronecker delta function, and 5° (a; — a;') the Dirac delta function. 
The three-dimensional PS of the universe is defined through 

(5(fc)5*(fe')) = P3D(fc)vO. (4) 

where (. . .) stands for an ensemble average. One may define an observed PS, p3D(fe) = \&{k)f jV , so that Pzuik) = (PsDik)). Note 
that a shot-noise term should be included if the PS is measured from discrete objects, and it is inversely proportional to the mean number 
density of the objects. We neglect the shot noise in this paper. For an isotropic universe, Psu{k) is a function of the length of k only, i.e. 
P3D(fc) = Paoik). The four-point function of S{k) is 

{5{k)t{k')S{k")S*{k"')) = l/2[P3D(fc)P3D(fc'0<5j^,„'5j^^„''' +f'3D(fc)P3D(fcO'5^,™''''5j^',™''] + V'r(fe,-fc',fe'',-fc'''), (5) 

where T is the trispectrum, and we have restricted the wavevectors to be in the same hemisphere so that the term P3D (fc)-P3D {k')5^ _„;/ S^, „„ 
does not appear. There is a redundancy in the variables of the tr ispectrum, wh ich has only six degrees of freedom arising from relative coor- 
dinates of the four points under the co nstraint of homog eneity 'Peebles"l98Cf). It is evident from the four-point function that the covariance 
of the three-dimensional PS is (see also lMeiksin & Whita>1999i : .Coorav & HtM2001i) 

a^{k, k') = {[Psoik) ~ P3D(fe)][p3D(fe') - P3D(fc')]> = Piuik)S^^^, + V-^T{k, -fc, k', -k'). (6) 

For Gaussian random fields (GRFs), the trispectrum vanishes, and a^{k,k') = {k)S^ „/ . In reality, the survey volume is always smaller 
than the observable univers e, and the survey geometry is more complex than the simple case we have assumed. These lead to a modification 
of the covariance for GRFs iFeldman. Kaiser & Peacockll994 



a\k,k')^Pio{k) 



Ju\x) e-'('=-'=''"da; 



J v?{x) Ax 



(7) 



where u{x) is a weight function that depends on the survey volume and geometry, and the shot noise is neglected. For the low-redshift 
cosmic density field, the non-Gaussianity can dramatically boost the elements of the covariance matrix through the trispectrum. 



3 ONE-DIMENSIONAL POWER SPECTRUM 

For a line-of-sight density that is along the x-3-axis and sampled at (xi, 2:2), the one-dimensional Fourier transform gives 

5{x±,m)^ J'' 5{x)e-^^'"^^^^^dx3^ ^ J2 <5(nx,n3)e2™^-^^/^, (8) 

*^ Tl ^ — — 00 

where the subscript _L signifies the first two components of a vector, i.e. x± — {xi,X2)- Similar to the three-dimensional PS, the one- 
dimensional PS is expected to follow 

{S{x^,m)S'{x^,n':,)) = PiDim)BS^^^„,^. (9) 

Substituting equation jsj in equation ^9} and making use of equation J4j, one finds the relation between the one-dimensional PS and the 
three-dimensional PS, 



© 0000 RAS, MNRAS 000, 000-000 



Covariance of One-Dimensional Power Spectrum 3 



PiTi{n-i) = ^ P3D(nx,n3), (10) 



Tl I — — OC 



which is a discrete analog of equation Q. 

Practically, one measures PS's of LOS densities sampled at some locations, e.g. Vi-d{x^, hj,) — \5{x_i_, k^)'^ /B. A simple estimator 
of the one-dimensional PS may be constructed by a spatial average over many LOS's, i.e. Vioik-i) = {Vi-d[x±, kz))x^, where (. . 
runs over all LOS's sampled in a single universe. To assess the performance of the estimator, two questions need to be addressed: (1) How 
does ViTiiki) relate itself to V3o{k)\ and (2) What is the covariance of Vi-oikz) with respect to Piuika)! The rest of this section answers 
the former, and Section|4|the latter. 

For simplicity, we assume that LOS's are sampled regularly in transverse directions at an interval of b = B /m, where m is an integer. 
Each LOS has a length of B, and x± = l±b with hjh — 0, . . . ,m — 1. The estimated one-dimensional PS is then 



m— 1 oc 2 ^ 

Ij_—0tIj_ — — C>O T1.j^,Jj^ — — oo 



^ S{n±,n3)6*(n±^ mjj_,n3), (11) 



where the equality, 

m — 1 

1- X ^ 27ri (n — n' )l/m cK • n i -i i /^ 

— ye =On,n'+mi J = 0, ±1 , . . . , ±00, (12) 



1=0 

has been used to obtain the last line in equation Jl H . It is easy to show with equation <11> that 

Pioika) = {Vioika)) = {ViD{x^,k3)). (13) 

Note that there is no spatial average on the far right side of equation J13> . It seems that only in the limit m ^ 00 do 'PiD(fc3) and Vsuik) 
follow the same relation as equations Q and JlOt . i.e. 'PiD(fc3) = X/JT --00 'P'-^o{k±, hi)/ B'^ , where we have assumed (5(co) = so that 
only ~ (Oi 0) terms contribute. In other words, to reduce the uncertainties of the recovered Vaoik) and P3D(fe), one has to increase the 
sampling rate in the transverse direction. 

Without losing generality, one may choose m to be even, so that equation <1 1> can be re-arranged into 

in/ 2 00 

Pin{n:i) = ^ ^ Vw{n±,na) = ^ ^ P3D(n_L, ns) + ^(na), (14) 

= — m/2 



where 



1 

p3D(n±,n3) = ^ 5{n±+mjj_,n3] 



(15) 



and 



m/2 



^(713) = -^ ^ ^ S{n^+mj^,m)S*[n^+m{j^ + j'^),n:i]. (16) 



—00 ■ii.±= — m/2 
' j', ^0 



Equation <I4t provides two equivalent views of the sampling effect. On one hand, the one-dimensional PS is a sum of the aliased three- 
dimensional PS, p3Y>{k), that is sampled by a grid of m x m LOS's. On the other hand, it is a complete sum of the underlying three- 
dimensional PS with an extra term ^(^3) that is determined by the sampling rate and properties of the density field. 

The discrete Fourier transform cannot distinguish a principal mode at |fc| ^ fcNyq from its aliases at fc ± 2fcNyq, k ± 4fcNyq, • ■ •, where 
feNyq is the sampling Nyquist wavenumber. For example, ^Nyq ~ it /b if one samples the field at an equal spacing of b. The alias modes 
will be added to the principal mode if they are not properly filtered out before sampling (for details, see lHocknev & Eastwoo(il98lh . Since 
the cosmic density field is not band-limited, aliasing can distort statistics of the field. The distortion on the three-dimensional PS cannot be 
quantified a priori, because it depends on relative phases between the principal mode and its aliases. The alias effect is less pronounced, 
if amplitudes of the alias modes are much smaller than that of the principal mode. Since the three-dimensional PS decreases toward small 
scales, a high sampling rate (or /cNyq) can suppress aliasing for modes with k <C /cNyq. 

If the Fourier modes of the density field are uncorrelated, the term A{ks) may be neglected even for a finite number of LOS's and, 
therefore, validate equation <I0> . Strictly speaking, A{k3) vanishes only as an ensemble average over many GRFs, but since there are so 
many independent modes in a shell of radius around k - especially at large wavenumbers ~ the summation in ^4(^3) will tend to vanish even 
for a single GRF. The cosmic density field is more Gaussian at higher redshift, so equation <10t mav be a goo d approximation then. At low 
redshift, however, the non-Gaussianity alters the statistics of one-dimensional fields (see e.g. I Amendolal 1 9941) . such that one might only be 
able to recover a heavily aliased three-dimensional PS from a sparse sample of LOS's. 
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4 COVARIANCE OF THE ONE-DIMENSIONAL POWER SPECTRUM 

The covariance of the PS is of interest because it tells us how likely an estimated PS is to represent the true PS and how much the modes on 
different scales are correlated. For the one-dimensional PS, the covariance is defined as 

fTm(fc3,fc3) = ([PlD(fc3) -PlD(fc3)][PlD(fc3) -PlD(fc3)]). (17) 

The covariance of the mean PS of N LOS's that are sampled in a single field can be expanded into a sum of pair-wise covariances between 
two LOS's separated by = x-'j_ — x''j_, e.g. 

1 ^ 

(^iDik:i,k':i) ^ afuik3,k'3;s^l), (18) 

i,i=i 

where 

aloik3,k's;s'l) = {[Viu{x{, k-,) ~ Piuik:i)]['Pii^ix'±, k'^) - Pm{k's)]) , 

and is the location of the jth LOS in the X1-X2 plane. 

For GRFs, the four-point function [equation ^5}] helps reduce the pair-wise covariance to 

1 

0-10(713, n3;s) = — ^ {S{n^,n3)S*{n'j_,ns)S{n'i,n':i)5*(n'",n's)) 



n I .Ti , ,71 .Tl. ' — — CX3 



xe 



27ri[(Ti._|_ — Ti'i )-a::'^ -\-{Ti'!—n'!')-x\]/B 



- Pir>{n3) Putin's) 



52 X! ■P3D(nx,n3)e^'''"-L"'/-^ S^^.n'^ = \^{s,n3)fS^,^„,^, (19) 



where the subscript and superscript are dropped for s. Fouri er transforms of ^(s,fc3) will give the three-dimensional PS and correlation 
function of the density field, and Pid (^3) — ^(0, fc3) (see also lviel et alj2002l) . Because of isotropy, the pair- wise covariance (Jin(fc3, ^3; s) 
and ^(s, fc3) depend only on the magnitude of the separation, i.e. (Jin(fc3, k'^; s) = aij^{ks, k'^; s) and ^(s, ^3) = ^{s, k^). 

If one LOS is sampled from each GRF, the variance of the measured one-dimensional PS Vio{k3) is ai-Q{k3, kg) = a^j-ilks, ^3; 0) = 
Pioiks), analogous to the three-dimensional case. If iV = LOS's are sampled in each GRF on a regular grid as in Section|3] e.g. s = 
{l± — l'±)b, the covariance becomes 



^ m— 1 ^ 00 m— 1 ^2 

iriu{ni,n's) = — ^ crlo[n3,n's-{l^-l'^)b]^ jjr^ ^ -P3D(n_L, n3)P3D (n^, 713) | ^ e^''''"^""^'''^''" 

I I .i'l — 71 I ,Tl' — — 00 I ^—0 



^ -P3D(nx,n3)P3D(nx +mj^,n3)(Sj^3_„^, (20) 

where we have used equation <12> to reach the last line. This result can be easily obtained using equations js) and <1 1> as well. It is seen that 
the summation in the last line of equation <20t runs over only one mode out of every A'^ modes in the Fourier space. If the three-dimensional 
PS were constant, the variance of the mean PS of the N LOS's would be A'^ times smaller than the variance of the PS of a single LOS. This is 
coincident with the theory of the variance of the mean. Since the cosmic density field is not a GRF in general, equation I20i is not expected 
to give an accurate estimate. 

If the A'^ LOS's are sampled randomly in each GRF, the summation over LOS's in the second line of equation ilOi should be re-cast to 

read 

00 N 

'^iDin3,n3) = j^^:^ 2^ P3D(n_L,n3)P3D(nx,n3) 2^ e ^-^ ^ ^ S„^ „,^. (21) 

n _\_,n'^— — oo j, 1 = 1 

Since x-'j_ is randomly distributed, the second sum in equation i2\\ tends to vanish for a large number of LOS's except when j = L Thus, 
one obtains 

cr?D(w3,74) - -^^^ X! -P3D(n±,n3)P3D(nY,n3)5^3,„/^ = — Pid(713)(5J^3,„^. (22) 

7T ,7T Y — ~ °0 

Again, it shows that the variance of PiD(fc3) is inversely proportional to the number of LOS's, but this is valid only for a large number of 
LOS's randomly sampled in a GRF. 

Through the Gaussian case one can find the lowest bound of the uncertainty for a measured one-dimensional PS and the minimum 
number of LOS's needed for a target precision. For example, to measure the one-dimensional PS of GRFs accurate to 5% on every mode. 
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one needs at least 400 LOS's. However, the cosmic density field has a far more complex covariance due to its non-vanishing trispectrum, and 
it can have a much larger variance of the measured one-dimensional PS. 

The trispectrum introduces an extra term to the covariance of the one-dimensional PS in addition to the Gaussian piece [equation <20H . 

i.e. 

o-iD(n3,ri3) = -gj ^ -P3D(nx,n3)P3D(nj_ +mj\,n3)5j^3_„/^ + -gTiD(n3,n3), (23) 

"i, Jj_ = -oo 

where 

1 ^ 

Tw{n3,n'^) = — ^ T{n,-n-mj,n',-n' -mf ), (24) 

n±, j j_, n'^, and j'^ are the transverse components of n, j, n', and j' , respectively, and js = j'^ — 0. In deriving equation <23> . we have 
made use of equations js} and <11> . 



5 LINE-OF-SIGHT LENGTH 

Observationally, the length of a LOS is always much less than the size of the observable universe. For example, in the case of a pencil-beam 
survey, the length is limited by the depth of the survey. For the Lya forest, it is often determined by the distance between the starting redshift 
of LyQ and Ly/3 lines. Damped Lya systems and other astrophysical or instrumental factors can break the spectrum into yet shorter chunks. 
Even if a very long LOS were obtained, one might still wish to measure the PS on smaller segments to avoid evolutionary effects. To account 
for the length of LOS's, the one-dimensional PS and its covariance must be re-formulated. 
For a LOS from (x±, 0) to {x±,L), its Fourier transform is 

5(a=x,n3)= / <5(xx,:r3)e-2-«^^3/^dx3 = -^ <5(n)e^— ^^^/^ " \ (25) 



n — — oo 



SO that 



(5(a;x,n3)(5*(a;x,n3)) = -g3 ^ P-iT,{n)w{n-i,n3)w{n3,n3), (26) 

Tl — — OO 

where 

t - \ / ixn3sin[(n3L/B -^3)71] 
w(n3,n3) = (-1) 7 • (27) 

The function w(n3, 723) is a window function that mixes Fourier modes of the density field along direction into the one-dimensional 
Fourier mode at fi-j,. Because 10(713,723) and 7.0(723,77.3) are not orthogonal to each other, the term {S{x^,hs,)5*(x^,n'-^)) is no longer 
diagonal with respect to 723 and 72,3. However, it remains diagonal-dominant. 
One may define the observed one-dimensional PS at a; x as 

PiD(a;x,n3) = (i5(xx,n3)|')/L. (28) 

The ensemble-averaged one-dimensional PS is 

Piv{na) = {Piuix_i_,h3)) = — ^ P3D(n)7;;^(7i3, 723). (29) 

n — — 00 

For GRFs, the covariance of Viuins) is 

o-iD(n3,"-3) = ([7'iD(a;x,?23) - -PiD(".3)][PiD(a;x,n3) - PiD(n3)]) 

= — ^ P3D(n)P3D(n')77;(n3, 723)10(723, 723)10(723, n3)77;(723, 223). (30) 



PS 



As expected, the covariance is not diagonal because of the window function 20(223, 723) 



6 Af-BODY TESTS 

Even at moderately high redshift, the cosmic density field is already quite non-Gaussian on scales below 10 ftr^Mpc. When projected in one 
dimension, the small-scale non-Gaussianity will obviously affect the measured one-dimensional PS on much larger scales. The non- vanishing 
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1 23401 23401234 

kg (7r/16 h Mpc-i) 

Figure 1. Normalized pair-wise covariance matrices C{kz, k'^; s) [see equation <32H in grey scale. The covariances are calculated over 2000 pairs of LOS's 
sampled at a fixed separation s from the simulation B512 at ^ = 3 (lower panels) and from 2000 GRFs (R512) that have the same box size and three- 
dimensional mass PS as the simulation (upper panels). The grey scale extends to —0.1, which corresponds to white. 



trispectrum introduces an extra term, riD(A:3, fcg), to the covari ance of the one-dime nsional PS [see equation <23H . Although it is possible to 
derive an approximate trispectrum based on the halo model (e.g. lCoorav & HiJ200ll) . the one-dimensional projection, unfortunately, obscures 
the contribution of the trispectrum to the one-dimensional PS. Therefore, numerical simulations are necessary for the study. 

Three A^-body simulations of 256'' cold dark matter (CDM) particles are used t o quantify the covariance of the one-dimensional PS. 
The model parameters are largely consistent with WMAP results JSpergel et all2003h . e.g. (fJ, fib, f^A, h, as, n) = (0.27, 0.04, 0.73, 0.71, 
0.85, 1), where Q is the cosmic matter density parameter, Qh the baryon density parameter, the energy density parameter associated 
with the cosmological constant, erg the rms density fluctuation within a radius of 8 h^^Mpc, and n is the power spectral index. The box 
sizes of the simulations are 128 h^^Mpc (labelled as B128), 256 h^^Mpc (B256), and 512 fe~^Mpc (B512). The baryon density param- 
eter f^b is used only for the purpose of calculating the transfer function using LINGER iMa & BertschingeJl 19951) . which is then read by 
GRAFIC2 (Bertschinger.,2001.) to g enerate the initial condition. The CDM particles are evolved from z = 44.5 to the present using GADGET 
iSpringel. Yoshida & Whitel200ll) . Additionally, GRFs are generated with the same box sizes as those of A''-body simulations for comparison 
and are labelled as R128, R256, and R512, respectively. 

The simulations produce snaps hots at z = 3 and 0. For each snapshot, the particles are assigned to a density grid of 512'' nodes using 
the triangular-shaped-cloud scheme iHocknev & Eastwoodll98ll) . LOS's are then sampled on the xi-X2 plane, and each LOS is extracted 
with full resolution along the xs-axis. Because the density grid has a finite resolution, both the coordinates and the wavenumbers are discrete 
and finite. The equations in Sections|2H5|need to be modified accordingly, so that they do not sum over non-existing modes. For example, 
equation <10> becomes 

M/2 

PMns) = ^ J2 PMn±,n3), (31) 

Ti^ = -A//2 

where M = 512 is the number of nodes of the density grid in transverse directions. 

For brevity and the purpose of comparing the covariance, we introduce the following three covariance matrices: (1) The normalized 
pair-wise covariance between two LOS's, 

C(fc3,4;s) = '^?D(fc3,4;s)[PiD(A:3)PiD(fc;5)]"\ (32) 

where s is the transverse separation between the two LOS's. (2) The normalized covariance of the estimated one-dimensional PS, 

C(fc3,4) ='^?D(fc3,fc3)[AD(fc3)PlD(fc3)]"\ (33) 

where (jfn(fc3, ^3) is the covariance of the mean PS of TV LOS's, e.g. equation jl 1> . and we have omitted its dependence on the number of 
LOS's. (3) The reduced covariance 
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Figure 2. Diagonal elements of the covariance matrices C'{k:j, k-j; s). The horizontal soUd Hne in the left panel marks the rms value of the off-diagonal 
elements in the s = 8 h^^Mpc case. Other lines are from direct summations of the three-dimensional mass PS using equation <19t . The symbols are 
measured from GRFs R512 (left) and the simulation B512 at z = 3 (light). For GRFs, C(fc3, k^; 0) is expected to be unity. 



C{k3, k'-,) = C{k3, fc3)[C(fc3, fc3)C(fc:^, k'^)]-^^^ (34) 

For GRFs, all these covariances are diagonal in matrix representation, if the length of LOS's is the same as the size of the simulation box. 
In addition, we have C(fc3, fc3; 0) — 1 and C(k3, fcs) = A'^"^, where is the number of LOS's that are sampled for estimating the one- 
dimensional PS. The advantage of (7(^3, fcg) is that Cik^, ks) = 1 for all fields, so that they can be compared with each other in a single 
(grey) scale. 

The normalized pair-wise covariance C(fc3,fc3;s) is shown in Fig.Qfor GRFs R512 and the simulation B512 at z = 3. The co- 
variances for GRFs are averaged over an ensemble of 2000 random realizations, while those for the simulation are averaged over 2000 
pairs of LOS's from a single field. The behaviour of the covariances is consistent with the expectation. Namely, C{kz, ^3; 0) ~ ^1 and 
C(fc3, fc3; s) decreases as the separation s increases. The simulation does deviate from GRFs because of the non-Gaussianity, which increases 
the variance C{k3, k^-^s). Fig.|2|compares the diagonal elements C{k3, k^; s) with those calculated using equation <19> . We note in passing 
that the expected values of C{k3, kz; 8 /j-^Mpc) are so close to that they are even below the rms value of the off-diagonal elements of 
C(k3, k'z;8 /i~^Mpc) for 2000 GRFs. Hence, it is practically difficult to recover three-dimensional statistics from 0-10(^3, fcs; s) or ^(s, ks) 
if the LOS's are too far apart. In theory, the modes of two LOS's are always correlated as k^ — > 0, regardless of their separation. However, 
large-scale LOS modes are greatly affected by small-scale three-dimensional perturbations, which are nearly-uncorrelated for large sepa- 
rations. As such, the correlation for s > 8 ft~^Mpc is so weak that it will not be easily detected against statistical uncertainties even for 
a sizable sample of 2000 pairs of LOS's. Hence, Figs.QandlJIsuggest that LOS's sampled in a single c osmic density fie l d are p ractically 
independent of each other as long as they are well separated. This condition is easily met. For example. iMcDonald et alJ )2004l) use 3035 
z > 2.3 quasar spectra from the Sloan Digital Sky Survey (SDSS, over 2627 square degrees) to determine the Lya flux PS. Assuming 
WMAP parameters, one will find that the mean separation between LOS's is ^ 60 comoving h^^Mpc such that the covariance arisen from 
the correlation between LOS's can be safely neglected. 

Fig.|3|shows the one-dimensional PS measured by averaging over 64 and 1024 LOS's from the B128 simulation. Although the mean 
PS of all the 512^ LOS's a grees with the result of a direct summation of the three-dimensional PS using equation <31> . the deviation of the 
PS for any particular group of 64 or 1024 LOS's is substantial, especially at 2 = 0. The variance is smaller at z = 3 than at 2: = because 
the cosmic density field is more Gaussian earlier on, and it is roughly inversely proportional to the number of LOS's. This can be seen better 
by comparing the lower right panel of Fig.|4|with that of Fig.|5| However, even at 2: = 3 the variance of the one-dimensional PS is still much 
higher than the variance for GRFs, 7V^^Pio(fc3), which indicates a heavy contribution from the trispectrum. The formulae in Sections|3Hl| 
often assume that L OS's are sampled on a grid with fixed spacing, which may not be applicable to realistic data such as inverted densities 
from the Lya forest iNusser & HaehnelJl999l:lzhaij2003l) . Therefore, we sample the LOS's in two ways in Fig.|3| grid sampling and random 
sampling. Since no significant difference is observed, random sampling can be safely applied in the rest of this paper. 

The normalized covariances C{k3, k's) and C{ks, k'^) are quantified in Figs.|4|(iV = 64) and|5l(iV = 1024). The left column in each 
figure is the covariances of the spatially averaged one-dimensional PS from a single GRF that has the same box size and three-dimensional 
PS as its corresponding simulation. Clearly, the covariances based on spatial average are nearly diagonal with unity diagonal elements. This is 
in agreement with the expectations based on ensemble average for GRFs and is consistent with the ergodicity argument. The middle column 
is similar to the left column except that the density fields are from simulations at 2 = 3. The modes in the simulated density fields are 
strongly correlated, so that the covariances are no longer diagonal. In other words, the trispectrum is non-vanishing for the cosmic density 
field, as it is the only term that contributes to off-diagonal elements in the covariance. The cosmic density field becomes so non-Gaussian at 
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Figure 3. Estimated one-dimensional mass PS's PidC^s)- The PS's are measured by averaging over 64 (liglit grey) and 1024 (dark grey) LOS's that are 
drawn from the B128 simulation. Shaded areas mark Icr dispersions of the PS's among 2000 (light grey) and 256 (dark grey, as 256 X 1024 exhausts all the 
512^ LOS's) distinct drawings. Circles are the mean PS's for dark grey areas, i.e. the mean of 512^ LOS's. Solid hues are results of a direct summation of the 
three-dimensional mass PS using equation 1311 . The LOS's are sampled on a grid with fixed spacing in the left panels but drawn randomly in the right panels. 

2 = that grey scale figures of the covariances will not be readable. Hence, we only plot four cross sections perpendicular the diagonal with 
Q = (fcs + fc3)/(7r/4 Mpc~^) for B 128 and B256 in the right column. The dominance of the diagonals suggested by these cross sections 
is actually weaker than that in the middle column, which can be seen by contrasting the cross sections for B 128 at z = in Fig.|4]with that 
forB128atz = 3inFig.|6l 

The non-Gaussianity is reflected not only in the correlations between different modes but also in the variance of the one-dimensional PS, 
i.e. the diagonal elements of the normalized covariance C(fe3, fca). The lower right panels of Figs.|4|and|5|compare C(fc3, fca) for five different 
density fields. The variance from the simulation B 128 is orders of magnitude higher than the Gaussian value of N^^ Plj^iks), and it grows as 
the non-Gaussianity becomes stronger toward z = 0. As a result, the sample variance error estimated for GRFs is much lower than what can 
be actually measured from the cosmic density field. According to equation <23> , both aliasing and the trispectrum contribute to the variance. 
Since the GRFs have the same three-dimensional PS and are sampled in the same way as the simulations, their near-Gaussian variances 
of spatially averaged one-dimensional PS's suggest that the contribution of the aliasing effect is negligible. Comparisons of C{hi, k-j,) for 
the same density field but with different sizes of sample {N = 64 and 1024) confirm the observation in Fig. |21 that the variance of the 
one-dimensional PS scales roughly as A'^"^. This means that even though the non-Gaussianity drives up the sample variance error of the 
measured one-dimensional PS, one can still reduce the error by sampling a large number of LOS's. 

The statistics of a large sample may sometimes be dominated by a small fraction of extreme cases. It is inevitable that some of the LOS's 
in our sample go through high-density halo structures, which may have significant contributions to the covariance matrices. To determine if 
a LOS falls in this category one has to assess the density profile intercepted by the LOS. We follow a much simpler approach by imposing 
a threshold to the mass column densities pcoi of the LOS's. The one-point distribution function of column density contrast pcoi/pcoi from 
the simulation B128 (upper left panel of Fig.|6| peaks at Pcoi/Pcoi = 1 as expected, and there clearly is a long-tail distribution (as compared 
to Gaussian or log-normal distributions) of Pcoi/Pcoi ^ a few. We re-calculate the covariance C{k-i, k'^) with a selection criterion that the 
column density of each LOS Pcoi/Pcoi < 3. This criterion excludes 318 LOS's from the selection of LOS's for the statistics. As a total of 
2000 X 64 LOS's are randomly selected from 512^ candidates, only half of the 318 LOS's, i.e. 0.1% of the sample, would be selected without 
the criterion. The result is shown in the upper middle panel of Fig.|6| which is visually similar to the B128 z — S panel in Fig.|4] However, 
the cross sections of C{k3, k'g) with the selection criterion (lower middle panel of Fig.|6} do demonstrate a several-percent reduction of the 
off-diagonal elements in contrast to those without the selection criterion (lower left panel). In other words, the 0.1% Pcoi/pcoi ^ 3 LOS's 
contribute at least 10 times as much to the covariance per LOS as pcoi/pcoi < 3 LOS's do. It is also useful to have a comparison of the 
normalized covariance C(fc3, k'^) because, unlike the reduced covariance Cik^, k'j) that only tells the relative strength between the diagonal 
and off-diagonal elements, C{ks, k'^) reflects absolute values of the covariance up to a factor of the product of PS's at fcs and k'^. Fig.0 
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Figure 4. Reduced covariances C{kz, k'^) [see equation 134H averaged over 2000 groups, each of which consists of 64 LOS's {N = 64). For each panel in 
the left, the LOS's are randomly drawn from a single GRF that has a box size of 128 h^^Mpc (R128), or 256 h^^Mpc (R256), or 512 h^^Mpc (R512). The 
GRFs have the same three-dimensional mass PS as their corresponding simulations at 2: = 0, but note that (7(^3, ^13) is independent of redshift for GRFs. 
Similarly, the middle column is for simulations at 2 = 3. The covariances are shown in a linear grey scale with black being 1.2 and white less than or equal 
to 0. At 2 = 0, the normalized covariances become much less diagonally dominant than the B128 z = 3 panel in the same grey scale, so only 4 cross sections 
along Q = (fcs + k'^)/{TT/4 h Mpc^^) are plotted for B128 and B256. Diagonal elements C(fc3, fcs) [see equation ^] of R128 (squares), B128 at 2 = 
(circles), B128 at 2 = 3 (diamonds), R512 (crosses), and B512 at 2 = 3 (pentagons) are shown with a multiplication factor in the lower right panel. For 
GRFs, NC{k:i,k^) = 1. 

shows the diagonal elements of the normalized covariance for B128. It is seen that the variance C{ks,ks) of the one-dimensional PS is 
indeed reduced by roughly a factor of 2 on all scales when the 0.1% Pcoi/Pcoi ^ 3 LOS's are excluded. 

For a fixed number of particles, the simulation box sets a cut-off scale, below which fluctuations cannot be properly represented in 
low-density regions such as the Lya forest. In other words, the number of particles and the size of the simulation determine the highest- 
wavenumber modes that are included in calculations of the one-dimensional PS and its covariance. Since the non-Gaussianity is stronger 
at smaller scales, a larger simulation box cuts off more small-scale fluctuations and may cause the correlation to appear weaker in Figs.|4| 
and|5| To test this, we assign the density field of the simulation B128 at 2 = 3 on a grid of 128'' nodes. The spatial resolution is the same 
as the simulation B512 on a grid of 512'' nodes. The covariance C{k3, k'j) is calculated in the same way as those in Fig.|4|but with fewer 
groups of LOS's. Each group still has 64 LOS's. The results are shown in the right column of Fig.|6|and in Fig.0 Evidently, there is a 
significant reduction of the correlations between different modes as well as the variance of the one-dimensional PS. Indeed, the variance 
from the low-resolution calculation is close to that from the B512 simulation. Thus, the apparent resemblance between B512 and the GRF 
R512 in Figs.|4|and|5|is mostly due to the low resolution of the large-box simulation. 
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It is expected from equation <30> tliat thie covariance matrix will not be diagonal if the length of LOS's is less than the size of the 
simulation box. Fig.|8|shows the covariances that are calculated in the same way as those in Fig.|4]for the GRF R512 and the simulation 
B512, except that each LOS is only 128 /i^^Mpc long. The effect of the length is not visible for the GRF, but it does increase the correlation 
between different modes and doubles the variance of the one-dimensional PS for the simulation B512. For real observations, the LOS length 
is always much less than the size of the observable universe, so that the window function w(k:i,, k-^) in the LOS direction will cause stronger 
mixing of modes and more pronounced increase of correlation and variance. 



7 CONCLUSIONS 

We have investigated the covariance of the one-dimensional mass PS for both GRF s and simulated density fields. The non-linear evo- 
lution and the non-Gaussianity of the cosmic den sity field on small scales (see also Izhan. .Tamkhedkar & FanslboOll ISmith et alJbOQgt 
IZaldarriaga. Scoccimarro & HuJ200.3l:lzharJ2003h have caused the correlation between the fluctuations on different scales and increased the 
cosmic variance of the one-dimensional PS. Because of this, large number of LOS's are needed to accurately measure the one-dimensional 
PS and recover the three-dimensional PS. 

The length of LOS's introduces a window function in the line-of-sight direction, which mixes neighboring Fourier modes in the cosmic 
density field. Fig.|8|has demonstrated for simulations that the length of LOS's (a quarter of the simulation box size) does affect the covariance 
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Figure 6. The effect of the non-Gaussianity and resolution on the covariance. The left column shows the distribution function of column density Pcol/Pcol 
(upper panel) and cross sections of the reduced covariance C{kz, fcg) along Q = (fcs + fc^)/(7r/4 /i Mpc"^) (lower panel) for B 128 at 2 = 3. The middle 
column shows C{k:j, fcg) (upper panel) and its cross sections (lower panel) for the same simulation output but with an exclusion of 318 LOS's that have 
Pcol/ Pcol 3- The reduced covariances in the first two columns are calculated in the same way as in Fig.El i.e. with 2000 groups and N = 64. Thus, about 
half of the 3 1 8 LOS's would be selected without the criterion p^^i / p^^i < 3. The right column is similar to the middle column, but the density is assigned on 
a grid of 128^ nodes. All the 128^ LOS's are selected and divided into 256 groups with A'^ = 64. The grey scale is the same as Fig. 141 
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Figure 7. The effect of the non-Gaussianity and resolution on the variance, i.e. diagonal elements of the covariance cr^p (fca, k'^). Diamonds are the normalized 
variance C{fc3 , fcs ) (multiplied by A'^) for B 1 28 at ^ = 3. Circles correspond to the middle column of Fig.|6] which imposes the selection criterion Pcol / Pcol < 
3. Squares are from the low-resolution calculation in the right column of Fig. 1^ 
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Figure 8. The same as the bottom row of Fig.|4| but with length of LOS's L = 128 h'^Mpc. 

of the one-dimensional PS. The covariance of the observed one-dimensional PS will receive more contributions from this effect, because in 
practice the length of LOS's suitable for study is always much less than the size of the observable universe. 

One may reduce the cosmic variance by binning independent Fourier modes. However, since the modes of the cosmic density field are 
strongly correlated, binning will be less effective in reducing the sample variance error. On the other hand, the non-Gaussian behaviour of 
the covariance provides important information of the field such as the trispectrum. 

Our analysis shows that even with 64 LOS's at z = 0, the cosmic variance on the P S can be several times o f the PS itself (see 
Fig.|4}. Thus, the extra power on the scale of 128 /i~^Mpc (SI = 1, SIa = 0) measured by ISroadhurst et all il99Ch from 4 samples of 
pencil-beam galaxy surveys can still be consistent with standard cosmogonies (with considerations of effects su ch as redshift distortion 
and non-Gaussianity, see e.g. lKaiser&Peacocldll99lL IPark & Gotllll99lLlvan de Wevgaerll99ll : lAmendolj|l994 . One may also infer the 
one-dimensional density field and measure the PS from the Lya forest. Since the Lya flux is the direct observable, many works have been 
based on fl ux sta t istics, which are then used along with simulations to constrain cosmology. Using A^-b ody simulations and the pseudo-hydro 
technique, IZhari y003|) find that the mapping between the nearly-Gaussian Lya flux PS IZhaj2004l) and the underlying one-dimensional 
mass PS (in redshift space) has a much larger scatter than that between the flux PS and the one-dimensional linear mass PS. This is expected 
because the one-dimensional PS of the cosmic density field has a variance much higher than the Gaussian variance. Therefore, one can 
measure the Lya flux PS precisely with a small sample of LOS's, but to recover the mass PS to the same precision more LOS's are needed. 
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